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FOREWORD 


For the past twenty-four years, the Naval Surface Weapons Center has been 
involved in a research and consulting effort for both the Energy Research and 
Development Administration and the Nuclear Regulatory Commission concerned 
with the study of reactor vessel response to hypothetical core accidents and 
other types of dynamic loading events. As part of this effort, analyses of jet 
forces produced by ruptured steam process pipes on neighboring walls of reactor 
buildings have been performed. This report presents the computed results for 
impact loads from the ruptured process pipe on the inner wall of a surrounding 
concentric guard pipe for geometries consistent with reactor plant design 
guide!ines. 
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CHAPTrR 1 
INTRODUCTION 


1.1 BACKGROUND 

Since 1956 the Naval Surface Weapons Center has been involved in a research 
and consulting effort for both the Energy Research and Development Administration 
and the Nuclear Regulatory Commission (NRC) concerned with the study of reactor 
vessel response to hypothetical core accidents and other types of dynamic 
loading events. In an earlier effort, NSWC provided to NRC preliminary calcula¬ 
tions of pressure loads produced on neighboring walls by steam pipe ruptures. 

The task reported herein is a continuation of these calculations for the 
purpose of determining the pressure loads very close-in to the site of steam 
pipe rupture. Specifically, a series of hydrocode calculations were done to 
determine the fluid flow/pressure fields in regions bounded by an inner steam 
pipe and an outer guard pipe following simulated ruptures of the high-pressure 
inner steam pipe. 

1.2 PROBLEM STATEMENT 

A process pipe, filled with saturated steam at 8.3 MPa (1200 psia), was 
assumed to rupture instantaneously in either the circumferential or longitudinal 
direction (See Figure 1). The steam flowed out from the rupture, driving 
ambient air ahead of it, and impacted the inside surface of a guard pipe which 
surrounded the process pipe. The desired result was the pressure distribution 
along the guard pipe. 

The analysis was performed in two phases. In the first phase, the proper 
equation of state for high-pressure saturated steam was evaluated, the appropriate 
finite difference grid scale for the problem geometry and flow conditions was 
determined, and a particular hydrocode for performing the calculations was 
selected. The second phase of the task determined the time-dependent internal 
pressure distribution on the, guard pipe as a function of (a) direction of the 
pipe rupture (circumferential or longitudinal), (b) width of the pipe rupture, 
and (c) separation distance or gap between the process and guard pipes. The 
parameters for the seventeen process/guard pipe configurations that were 
computed are listed in Table i. 
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FIGURE 1. GEOMETRY OF BREAKS 
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CHAPTER 2 

ANALYTICAL SOLUTIONS 


For a free compressible jet expanding to ambient pressure, the total jet 
thrust for an ideal gas as given by Reference 1 (slightly modified) is 


T/A„ - C.P A - P 
e 10 * 

where T » Jet thrust = Total force on a large normally-impacted plate 
A e » Effective area of process pipe break 
Pq * Pipe reservoir pressure 
P^ *= Ambient pressure 
C, - (1 + y)(2/(y+ 1)) y/(y ' 1) 
y - Specific heat ratio 

The coefficient C, is a weak function of v» being 1.255 for y ~ 1.3 and 1.229 
for y - 1.1058. 1 

This solution, or any other such simple solution, is not adequate here 
because 

(1) The gap between the pipes is relatively narrow in comparison with the 
pipe dimensions, and backpressure on the outside of the process pipe must be 
accounted for in a force balance; 

(2) The non-uniform filling of the gap between the process and guard 
pipes with steam and compressed air makes ill-defined; 

(3) The nozzle discharge coefficient is not known a priori to obtain the 
effective value of A g for the break in the process pipe; 

(4) Time-dependent pressures must be considered here initially to obtain 
the early-on transient airshock/steam shock contributions to the guard pipe 
loading; 

*Moody, F, J., "Prediction of Blowdown Thrust and Jet Forces," ASME Paper 69-HT-31, 
ASME-AIChE Heat Transfer Conference, Minneapolis, Minn., Aug 1969. 
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(5) The distribution of pressure along the guard pipe inside surface, 
not just the total force, is needed for future structural response calculations. 

The total force calculations are useful, however, for comparing with 
integrated results obtained from the hydrocode computations. These comparisons 
are discussed later. 


V 







NSWC TR 80-229 


CHAPTER 3 

HYDROCODE SOLUTIONS 


This was a problem well suited for an Eulerian (rather than Lagrangean) 
hydrocode since there is much distortion in the flow. In the initial phase 
of the analysis two hydrocodes were used, TUULI and CSQ. 

TUULI* is a two-dimensional Eulerian hydrocode, written in 1975 at NSWC 
by D. Lehto. Full documentation is not yet available. The code is based on 
the fluid-in-cell (FLIC) method (REFERENCE 2) where the calculation is done 
in four steps per time cycle: 

1) The accelerations are calculated from the pressure gradients and new 
provisional velocities are calculated (without convection). 

2) The provisional velocities are used to calculate the pdV work done on 
each zone; this gives provisional internal energies (still without convection). 

3) The material transport (convection) is done with the provisional 
velocities and energies. 

4) Any adjustments ..jeded to conserve both energy and momentum are made. 
Any kinetic energy correction needed for momentum conservation is taken from 
(or added to) the internal energy. These adjustments are necessary because 
the flow mixes dissimilar flows from adjacent zones, and both kinetic energy 
and momentum cannot be conserved in a mixing process. This fourth step is 
explicitly done in TUULI because internal energy is transported; it is 
implicitly done in the original FLIC code because total energy is transported. 
This is a trivial arbitrary choice. 

TUULI handles two materials (here, steam and air) and has an option for 
assigned inflow in any chosen zones (used extensively in these calculations). 

The grid is composed of fixed rectangular zones. Shocks are handled by the 
quadratic artificial viscosity method. 

-- 

Gentry, R. A,, Martin, R. E., and Daly, E. J., "An Eulerian Differencing 
Method for Unsteady Compressible Flow Problems," J. Comput. Phys . 1_, pp 87-118, 
1966. 

♦Until now, this code was called TUTTI, However, this is also the name of a Los 
Alamos Scientific Laboratory equation-of-state program (J. Appl. Phys. 5^ (10) 
5360 (1980)). Thus, the name change: TUTTI to TUULI. 
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CSQ is a two-dimensional Eulerian hydrocode written at the Sandia Laboratories. 
Documentation for this code is provided in Reference 3. In this method, the 
calculation is performed in two main phases: 

1) The finite difference analogs of the complete two-dimensional Lagrangean 
equations are solved during each time cycle. 

2) At the end of the time cycle, the code rezones the mesh back to the 
original configuration. 

The net result of the rezoning is an Eulerian calculation. The version of 
CSQ described in Reference 3 handles only two materials in any chosen zone. 

The grid is composed of fixed rectangular zones. Shocks are handled by the 
quadratic artificial viscosity method. 

The computations were performed on a CYBER 176 at the Air Force Weapons 
Laboratory in Albuquerque, New Mexico over telephone lines from NSWC. 

3.1 ZONE SIZE SELECTION 

Both TUULI and CSQ face the same limitations on available computer size 

and cost of computer time. Care was taken to choose a calculation mesh just 

fine enough to give the desired accuracy; this was done by simply trying 

progressively finer meshes until the overpressure loading on the guard 

pipe inside wall no longer changed significantly. The final mesh selected 

was 0.315 zones/mm (8 zones/in) within the channel formed by the break after 

it was determined that reduction of the mesh size down to 0.630 zones/mm 

(16 zones/in) produced no significant effect on the calculated flow. 

The mesh was nonuniform and varied for each problem geometry, being fine in 

the channel and near the impact area on the guard pipe and progressively 

coarser with increasing distance from the region of interest. 

» 

3.2 ZONE SHAPE LIMITATION 

Both hydrocodes are limited to rectangular zones. This makes the 
longitudinal-break problem awkward, because the concentric circles needed 
to represent the problem cross section would have to be made up of rectangular 
steps (See Figure 2). For the longitudinal-break calculations, the pipes were 
straightened out into parallel planes as shown in the figure; this is expected 
to be a good approximation since the separation distances between pipes are 
relatively small compared with the pipe radii for the geometries investigated 
here. No problem arose in setting up the mesh with rectangular zones for the 
circumferential-break calculations because all bounding surfaces could be 
represented by straight lines. The computing mesh for the flat-plate 
approximation for the longitudinal break shown in Figure 2 looks similar to 
the mesh for the circumferential break; however, the latter has axial symmetry. 


3 

Thompson, S. L., "CSQ -- A Two Dimensional Hydrodynamic Program with Energy 
Flow and Material Strength," Sandia Labs. SAND 74-0122, Aug 1975. 
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FLAT-PLATE APPROXIMATION 
FIGURE 2. GEOMETRY FOR LONGITUDINAL BREAK 
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3.3 EQUATIONS OF STATE 

When the steam pressure drops as it leaves the ruptured process pipe, 
condensation takes place and the resulting two-phase flow expands isentropically. 
For example, the isentropic exponent for equilibrium flow (i.e., equilibrium 
between the vapor and droplet phases) is 1.1058 (Table 2) for saturated steam at 

8.3 MPa (1200 psia).* However, it takes a finite time for the vapor to condense 
into droplets. If this relaxation were a relatively slow process with respect 
to the expansion, the expanding steam would continue to behave as a pure vapor. 

On the vapor side of the saturation line, the adiabatic exponent is 1,2592; 

a frozen flow would be expected tn maintain this exponent beyond the saturation 
line. The actual flow would lie between these two states during the expansion 
process. Referring to Figure 3 (taken from Reference 5), the Tl-curve is 
the equilibrium flow isentrope and the I'-curve segment below the saturated 
vapor curve (SV) is the frozen flow isentrope. The actual curve, nonequilibrium 
flow designated ACT, for the two-phase steam mixture which takes into account 
the kinetics of the relaxation/condensatiori process falls somewhere in between 
the TI and the I 1 curves as shown in the figure. 

The effect of nonequilibrium flow (delay in condensation) was bounded 
by performing hydrocode calculations with adiabatic exponents {y) for both 
the frozen and the equilibrium flow conditions. An ideal-gas equation of 
state. 


P • (y - 1)pE , 

was used to give the isentrope, 

Pp" Y = constant, 

where P is pressure, p is density, and E is internal energy. The results, 
given in Figure 4 for a typical circumferential-break geometry, indicate 
that nonequilibrium flow corrections to the pressure loading on the guard 
pipe opposite the pipe break are negligible. Equilibrium flow was assumed 
for the rest of the hydrocode calculations. 


♦Steam table (Reference 4) pressure-volume data were fitted (on a log-log 
scale) with straight lines to get the adiabatic exponents given in Table 2. 

^Keenan, 0. H. and Keyes, F. G., Thermodynamic Properties of Steam. Including 
Data for the Liquid and Solid Phase s, 1st'ed., Wiley,'New York, 1936. 

5 

Zel'dovich, Ya. B., and Raizer, Yu, P., Physics o f Shock Waves and High 
Temperature Hydrodynamic Phenomena, Volume II, ed. by U. b. Hayes and 
R. F. Probstein, Academic Press, New York, 1967. 
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TABLE 2 

saturated steam uata 


METRIC UNITS- 


PMESSURE 

(MPA) 

b.b*S 

b • 2 74 

b.blb 

SPECIFIC VOLUME 

(M3/MG) 

27.H2 

22.6V 

21.54 

i)EnSI fY 

(MG/M 3) 

O.OJb*b 

0.04426 

0.04644 

TEMPEkATURE 

(K) 

5S/.9B 

S70.S1 

573.39 

tNTROPY 

(J/KG.K) 

3blb. 

S722. 

5699. 

ADIABATIC CONSTANT 

ON VAPOR SIDE 


1.2623 

1.25*2 

1.2777 

ADIABATIC CONSTANT 
ON 2-PHASE SIDE 


1 • 0866 

1.1056 

i .104b 


t NGL I SH UNITS — 
KMESSURE 

(PS I A) 

1000 

1200 

1230 

SPEC IKIC VOLUMt 

(FT 3/LH) 

0.44bb 

0.3b IV 

0.3460 

Tt^PtHATUPE 

(F > 

344.67 

S67.22 

5 7 2 • 4 2 

ENTROPY 

(BTU/F.LB) 

1.3H*7 

1.3667 

1.3bl2 


TABLE 3 

tAlT CONDITIONS AT SLIT 


b.274 MPA <1200 PblA) SATURATED 5 T t AM , G AMM A = 1.1054 


bkeak wioth« 

b0.« MM (? 

IN). PIPE 

thickn tbs 

=50.b MM 

(2 IN) . 

distance FROM 


pakticle 

5UOND 

MASS 

MluPLANE OF 

PKtSSURE 

UtNsiTY 

VELOCITY 

sPtED 

FLUX 

BhEAK (MM) 

(MPA) 

(G/MM3) 

(M/S) 

(M/S) 

(G/M2/US) 

1.6 

2.613 

lb, 52 

6UV.H 

431.5 

*.464 

4.8 

2.616 

15.32 

hOb.5 

431.7 

* .444 

/.9 

2.61b 

15.32 

605. V 

431 . V 

*.404 

11.1 

2.617 

15.52 

601.1 

431 .B 

*.32* 

14.3 

2.606 

lb.52 

3*1.7 

430. V 

*. lbj 

17. b 

2.S73 

15.52 

573.2 

426.2 

b.b* 6 

20.6 

2.506 

15.52 

335,H 

422.6 

a. 316 

23.B 

2.403 

15.52 

463.7 

4l3.b 

7.19 7 

SONIC 


27.6b 

472.B 

*♦72.6 

12.05 
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FIGURE 4. EFFECT OF THERMODYNAMIC NONEQUILIBRIUM ON PRESSURE ON GUARD PIPE 
OPPOSITE THE BREAK 
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The air in the gap between the process and guard pipes was treated as 
an ideal gas with y = 1.4. The air in the gap was initially set to an average 
temperature of 422 K (300° F) for all the hydrocode calculations. This value 
was determined by performing a steady one-dimensional heat conduction analysis 
(see Reference 6, p. 37) for a steam pipe surrounded by an air gap, a guard 
pipe, and an ambient atmosphere with free convection. 

3.4 HYDROCODE SELECTION 


TUULI and CSQ give essentially the same results, as expected. This is 
demonstrated in Figure 5 which presents the pressure versus time history from 
each code at the guard pipe surface directly opposite the circumferential break 
for Run No. 4 listed in Table 1. The small oscillations are unresolved airshock 
reflections between the nose of the steam jet and the guard pipe wall. They 
are considered insignificant in terms of pipe response to the pressure loads. 


TUULI was the hydrocode 
because the authors are more 


chosen for the remainder of the calculations 
familiar with it than with CSQ. 


^Kreith, F., Principles of Heat. Transfer , 2nd ed., International Textbook Co., 
Scranton, Pa., 1965. 
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TUULI 





8.274 MPa (1200 PSIA) SAT. STEAM 
50.8 mm (2 in) CIRCUMFERENTIAL BREAK 
PROCESS PIPE O.D. - 0.914 m (36 in) 
GUARD PIPE I.D. - 1.016 m (40 in) 


TUULI RUN T7887 
CSQ RUN 3 


TIME(mt) 


FIGURE 5. TUULI AND CSQ CODE PRESSURE VS TIME ON GUARD PIPE OPPOSITE THE BREAK 
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CHAPTER 4 
HYDROCODE RESULTS 


The expected flow phenomena following an instantaneous break in the process 
pipe (for the circumferential-break geometry) are shown in Figure 6. A 
rarefaction wave recedes into the process pipe steam reservoir. A steam plume 
jets from the break in the process pipe into the air gap between the process 
and guard pipes. Sonic conditions are attained within the break and the flow 
becomes choked. An airshock precedes the steam jet and a stationary shock 
front (in the steam) is formed which stands off at some distance from the 
guard pipe inner wal1. 

In section 4.1 several elements of the flow phenomena -- rarefaction wave, 
sonic flow in the process pipe break, and discharge coefficient for the break -- 
are discussed. These elements allow a substantial simplification of the flow 
field calculations. The flow internal to the process pipe need only be computed 
once (with a hydrocode) for each initial steam flow condition. Following this, 
truncated flow field computations are then performed for the different process/ 
guard pipe geometries. The truncated flow field calculation involves specifying 
the steady outflow from the process pipe exit and then computing the flow field 
for the intervening space between the process pipe and the guard pipe. 

Section 4.2 discusses the following aspects of the flow field for the 
region between the process pipe and the guard pipe: 

1) The pressure loads on the guard pipe wall produced by the airshock 
preceding the steam plume; 

2) The transient initial pressure loads produced by the steam plume on 
the guard pipe wal1; 

3) The stationary shock in the vicinity of the guard pipe wall; and 

4) The maximum pressure loads on the guard pipe wall that are established 
by the steady flow of the steam exiting from the process pipe break. 

4.1 OUTFLOW FROM PROCESS PIPE 

RAREFACTION WAVE . A rarefaction wave runs from the break into the steam in 
the process pipe. For a circumferential break, this wave reflects from the pipe 
axis and part of it comes back into the break. Figure 7 shows the guard pipe 
wall reflected pressure versus time for such a calculation carried out far 
enough in time to follow the propagation of the reflected rarefaction wave to 
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FIGURE 7. EFFECT OF AXIAL REFLECTION OF RAREFACTION WAVE IN PROCESS PIPE 
ON GUARD PIPE PRESSURE 
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and into the break. The net effect on the guard pipe reflected pressure loading 
is a negligible pressure oscillation that arrives after an essentially steady 
flow pattern at the guard pipe wall has been attained. The shape of the 
pressure oscillation indicated in Figure 7 represents the interaction of the 
rarefaction wave with the stationary shock and the guard pipe wall -- the lack 
of definition of this pressure oscillation is caused by the coarse zone size 
used in this computation. 

SONIC FLOW IN PROCESS PIPE BREAK . Within the break, sonic flow conditions 
occur at the minimum flow cross section if the downstream pressure, p^, is 



where p u is the upstream (reservoir) pressure and y is the specific heat ratio 
of the gas. 

For a free jet, the downstream pressure is equal to the ambient pressure. 

For a confined jet such as discussed here, the downstream pressure is not 
defined a priori because of the pressure buildup in the confined space. This 
does not present a problem for the hydrocode calculation because it automatically 
simulates the correct flow, whether it is sonic or subsonic. However, if 
preliminary hydrocode calculations indicate that the flow in the break remains 
sonic for the duration under investigation (which turns out to be the situation 
here), then considerable savings in computer costs can be made by assigning 
the proper steady flow conditions at the break exit and calculating only the 
flow downstream of the break for the various pipe spacings. The downstream 
flow cannot send any signals through the sonic f1ov/ into the process pipe. 

DISCHARGE COEFFICIENT FOR PROCESS P I PE BREAK . The process pipe break is 
expected to have sharp edges. For"such" a nozzle," the flow contracts and 
does not fill the "nozzle" cross section (Figure 8). This effect is the 
familiar vena contracta of incompressible flow and is known tc occur for 
compressible flow as w ell (References 7 - 11). The area contraction 

^Bean, H. S., Buckingham, E. and M urphy, P. S., "Discharge Coefficients of 
Square-Edged Orifices for Measuring the Flow of Air," N.B.S.J. Res. 2, 
pp 561-568, 1969. 

^Stanton, T. E., "On the Flow of Gases at High Speeds," Proc. Roy. Soc . Ill , 
pp 306-339, 1926. 

^Perry, J. A., Or., "Critical Flow Through Sharp-Edged Orifices," Trans. Am. 

Soc. Hech. Engrs . 71 , pp 757-764, 1949. 

^Grace, H. P. and Lapple, C. E., "Discharge Coefficients of Small-Piameter 
Orifices and Flow Nozzles," Trans. ASME , pp 639-647, Jul 1951. 

^Arnberg, B. T., "Review of Critical Flowmeters for Gas Flow Measurements," 

J. Basic Engq ., pp 447-46'', .tec 1362. 
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FIGURE 8. ROUNDED VS SHARP-EDGED INLET 
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coefficient is typically between 0.6 and 1.0. This contraction allows the 
downstream pressure to have some influence on the outflow even though the flow 
is sonic, because signals can reach the sonic region by bypassing the sonic 
core of the flow. Because no data were available on the flow of steam through 
sharp-edged channels, the "nozzle" or process pipe break flow was determined 
as part of the problem. The hydrocode calculations gave a discharge 
coefficient of 0.77, compared with 1.00 expected for a rounded-inlet orifice. 
This discharge coefficient is used in calculating the outflow from the process 
pipe break exit for some of the truncated flow field calculations. 

The steam jet reaches the end of the 50.8 mm (2.00 in) long break (process 
pipe wall thickness) before the flow has expanded to fill the entire cross- 
sectional area for the 50.8 mm break width. This result is evident in Table 3 
which gives the steady-state outflow parameters, computed using CSQ, for the 
8.3 MPa (1200 psia) process pipe reservoir conditions. The mass flux near the 
centerline is 30% greater than the mass flux near the "nozzle" wall. The 
effect is also indicated schematically in Figure 8 for the "nozzle" with the 
sharp-edged inlet -- the figure indicates that the supersonic flow has not 
expanded to encompass the entire area of the "nozzle," Results such as given 
in Table 3 will be used to avoid calculation in every computer run of the 
interior of the process pipe and the interior of the break. 

TRUNCATED FLOW FIELD CALCULATIONS . To reduce computer time costs, the 
runs with 8.3 MPa (1200 psia) process pipe reservoir conditions (Table 1, 

Runs No. 1 - 12, 16, 17), were computed using a truncated flow field. The 
calculation of the flow in the interior cf the process pipe (the reservoir) 
and the interior of the break (the "nozzle") was replaced by the steady-state 
flow conditions of Table 3 at the break outlet since the effect of the rare¬ 
faction wave moving into the reservoir was shown to be negligible in Figure 7. 

For the 12.7 mm (0.5 in) wide breaks in these runs, the flow filled the 
break exit, so the exit conditions were obtained from analytical sonic flow 
calculations over 0.77 (the discharge coefficient determined for these sharp- 
edged breaks or "nozzles") of the break area expanded to the full exit area 
with the usual shock-tube area relations. The sonic conditions given in the 
last line of Table 3 apply across the entire exit plane of the 12.7 mm (0.5 in) 
wide breaks. 

4.2 PRESSURES ON GUARD PIPE WALL 

AIRSHOCK FORMATION BETWEEN PIPES . The exiting steam jet drives the air 
between the pipes ahead of it, and forms at first a bow airshock at the nose 
of the steam plume, and later an airshock that propagates ahead of the steam 
jet flow in the space between the pipes. 

Figure 9 presents a typical pressure distribution along the guard pipe 
inside wall showing the presence of the airshock. The airshock shown in 
the figure has an overpressure of 0.5 to 1.0 MPa (5 to 10 bars) and is not 
regarded as structurally significant; it is therefore handled with coarse 
zoning. 
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FIGURE 9. PRESSURE DISTRIBUTION ALONG GUARD PIPE, INCLUDING AIR SHOCK 
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In the remainder of this report, the discussion concentrates on the sharp 
pressure maximum on the guard pipe inside surface facing the site of the process 
pipe break. 

PRESSURE TRANSIENTS ON THE GUARD PIPE HALL . Before discussing the main results, 
the pressure transients due to shock reflections between the oncoming steam jet 
and the guard pipe are examined briefly. These are poorly resolved in the 
two-dimensional calculations (e.g., see the oscillations in Figure 5). To 
resolve these shocks, the zoning scale would have to be about an order of 
magnitude finer and the computer cost for a two-dimensional calculation would 
be quite high. However, an upper bound can be obtained with a finely-zoned 
one-dimensional (WUNDY code; Reference 12) Lagrangean calculation. The 
simple linear shock tube geometry and the relationship with the geometry of 
the real problem are shown in Figure 10. The resulting pressure-time record 
at the guard pipe is shown in Figure 11. The reflected shocks are not 
significant because: 

1) They are very short (about 10 ps) and carry litt'e impulse; 

2) They would be reduced by the two-dimensional expansion for the real 
problem geometry; 

3) They would be very much reduced in amplitude by the finite opening 
time of the process pipe break which is assumed to be instantaneous for these 
calculations. 

The opening time of a circumferential 50.8 nm (2.00 in) break in an 
infinite pipe with the parameters of Runs 1 through 12 is about 27 ms (Appendix 
A); the opening time of a longitudinal break is about 8 ms, calculated using 
a clamshell-type opening model with hinging at the side opposite the break. 

These long opening times mean that the steam jet that first reaches the 
guard pipe originated from a small slit and would be greatly attenuated in 
pressure because the distance to the guard pipe would represent many slit 
widths. 

These insignificant transient effects were ignored. The computational 
effort was directed towards determining the steady-state pressure distribution 
on tne guard pipe. The slit (process pipe break exit) was assumed to open 
instantaneously, which allowed rapid attainment of the same steady flow that 
would have been eventually reached with a slowly-opening slit. 

The steady-flow pressures reported in the next section represent over¬ 
estimates if the pipes are so short that the supply pressure drops before 
steady flow is attained. 
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lehto, D. and Lutzky, II., "One-Dimensional Hydrodynamic Code for Nuclear- 
Explosion Calculations," Naval Ordnance Laboratory, NCLTR 62-168, Mar 1965. 
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FIGURE 10». GEOMETRY OF REAL PROBLEM 


PLANE 

SYMMETRY 



FIGURE 10b. GEOMETRY OF LINEAR SHOCK TU3E PROBLEM 
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STEADY-FLOW PRESSURES ON GUARD PIPE WALL . The main results of this study 
are presented In Figures 12 - lf>. These figures present the steady-flow 
pressure distribution on the guard pipe inside surface in the vicinity of the 
process pipe break. The full distribution of computed pressure data including 
the airshock effects are listed in corresponding Tables 4-6. The longitudinal 
break maximum pressures are greater than the corresponding circumferential-break 
maximum pressures, because the latter flow has more geometric spreading (compare 
Figures 12 and 13 and see Figure 14). 

In both process pipe break configurations, the maximum pressure drops off 
rapidly with increasing process/guard pipe separation distances. 

A check was made on how the total force on the guard pipe plus the process 
pipe agrees with the thrust equation for a jet of outlet area A hitting an 
unconfined plate (Reference 1). Repeating the expression given in Chapter 2, 

T/A e - Cl P 0 - P» 

with « 1.229 (for y » 1.1056), P 0 = 8.274 MPa and P* = 0.101 MPa 
gives T/A fi - 1.017 x 10 8 dyne/cm 2 * 10.17 M/m 2 . 

This can be compared with the net outward force on both pipes for the 
parallel-plate model of the longitudinal break: 


F out‘/V fl - /V dfl 

gp Pp 

For small plate separations, this integration removes most of the effect of 
the shock running between the plates, because the shock pressure distributions 
are nearly identical on the two surfaces. For large plate separations, the 
shock contribution does not cancel. 

Figure 16 shows that the calculated net forces are close to those of the 
thrust equation. These integrals should not be regarded as precise, because 
they depend strongly on results in the rather coarsely zoned airshock region. 
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FIGURE 12. STEADY-FLOW PRESSURE ALONG GUARD PIPE FOR 50.8 mm (2 in) 
CIRCUMFERENTIAL BREAK 
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FIGURE 13. STEADY-FLOW PRESSURE ALONG GUARD PIPE FOR 60.8 mm (2 in) 
LONGITUDINAL BREAK 
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TABLE 4 

STEADY FLOW PRESSURES ALONG INSIDE Of GUARD PIPE DUE TO STEAM JET 
FRO M 50.8 MM (2 IN) RREAK IN PROCESS PIPE 


0,^14 M (36 IN) U.t>. PROCESS PIPE 
B.274 mRA <1200 wsia) SATURATtD STEAM 

Z = UISTmNCE ALUNO GUAM!) PIPE f HUH CENlTfcR OF bRtAK (MM). 
F=AbSOuUTE PRESSURE On bUAPU HlPt (MPA). 

SRaCE= 3PACE HEfRtEN RlPEs (mm). 


RKOd.NO. 
SPACE 
/. 

1.6 
4 ♦ H 
6,0 
11.2 

14.3 

17.5 
2U.7 
23.9 
27.0 
SU.3 
34.0 
38 • 0 

46.3 

47.3 

32.5 
3 i .5 
b2.5 

67.5 

72.3 
7 6.5 

6b,0 
43.0 
10b.0 
U9.0 
137.0 
136.0 

155.0 
220.0 
263,0 
323.0 
415.0 
523,0 
635.0 
743.0 
660.0 
980.0 
1100.0 
1230.0 
1370.0 
1320.0 


C IRCUMf- tP£N U AL RWtMK 


787 0 

7 a 71 

7672 

50.6 101.6 152.4 

H 

R 

R 

6 . Ob 

3.14 

1 .82 

6.03 

3.14 

1.60 

3.42 

3.12 

1.60 

3.7 7 

3,10 

1.61 

5.3 7 

3.06 

1.80 

3.31 

3.00 

1.74 

3.0 3 

2.*3 

1.77 

4.72 

2.6 7 

1.73 

4.34 

2.60 

1.74 

4.04 

2.71 

1.71 

3.6b 

2.62 

1.64 

3. JO 

2.51 

1.64 

2.42 

2 • 34 

1 .61 

2.3a 

«v 

• 

rv 

1 .5b 

2.23 

2.11 

1.50 

1,94 

1.46 

1.43 

1.70 

1 .63 

1,34 

1.48 

1.64 

1,33 

1.30 

l .53 

1.2 7 

1,11 

1,40 

1.21 

.43 1 

1.23 

1.14 

.313 

1.04 

1.0 7 

. 731 

.431 

• 463 

.724 

.734 

. 684 

.437 

.6 11 

.772 

.444 

.474 

.653 

,434 

.37 2 

,5 33 

.■•■*0 2 

. 304 

.421 

• h23 

.3 32 

.336 

.771 

.416 

.306 

. 733 

.454 

.38 0 

.7 14 

.443 

,434 

.7 00 

.311 

.406 

.686 

.437 

.292 

. -4b 

. 2b4 

. 163 

.333 

.14 7 

.112 

.18 7 

.10 7 

.102 

.111 

.10 1 

.10 1 

.102 

.101 

.101 

.10 1 

.101 

.101 


U HS1) 

% 


ONUlTUUlNAL 

BPt AK 

7684 

7866 

767 3 

12.7 

23.4 

30.3 

P 

P 

P 

8.56 

6.19 

6. 14 

9.^6 

8.16 

6. 13 

6.27 

6 » OS 

6,20 

8. U 

7.90 

6.03 

7.47 

7.67 

d.bb 

7. 72 

7.57 

3.39 

7.57 

7.01 

O.30 

b • 44 

b.5B 

4.97 

6.29 

6.13 

4.63 

3.67 

3.64 

4.28 

3.12 

3.15 

3.40 

4.71 

4.66 

3.31 

4.44 

4.20 

3.10 


7874 )uf3 


4,2 rt 
4.1H 
4.10 
4.0 1 
3.*2 
J.62 
3. 72 
3.Cl 

3.44 


3 c 
22 
05 
HH 
66 
2.46 
22 
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2 . 
1 . 


.71 

.30 

• 33 
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• 1 o 
.06 
9U0 
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1 74 
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3.62 

3.32 
3.31 
3. 16 
3.0 7 
3.00 
2.44 
2.66 
2.60 
2.64 

2.36 
2.42 
2.26 
2.10 
1 .41 
1.72 

1.32 

1.32 
1 . 1H 
1.04 
1.04 
.464 
. 740 
.474 
« 166 
.104 
.101 


2.7 2 

2.37 
2.07 
i.bl 
1.36 
1.34 
1.16 
l.Ul 
.663 
.771 
,7 63 
.6/6 
,4^4 
1.01 
.444 
.664 
.602 
.760 
.734 
.7 16 
.7 03 

• b21 
.416 
.204 
.117 
.102 
.lul 


3.31 
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2 » 03 
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1. 74 
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.664 
.666 
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.332 
.342 
.4 34 
.462 
.510 
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.111 
.102 
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. \ ik i 


2. 16 
2 . 14 

2.14 
2 . 16 
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2. 13 
2.10 
6.06 
2.03 
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1.47 
1 .43 
1 • n 7 
1.60 
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l . 66 
1.36 
1 » 3 i 
1.43 
1.3<* 
1.25 

1.1 J 
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• b 12 
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.374 
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.446 
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TAHLt b 

STEADY FLOW PRESSURES ALONG INSIDE OF GUARD PIPE DOE TO STEAM JET 
FROM 12.7 MM (0.5 IN) BREAK 11! PROCESS’PIPE 


Q.V14 8 (36 IN) 0.0. FWoCtSS PIPE 
8.274 MPA < 1200 PSIA) SATURATED STEAM 

Z = 01 ST ANCE ALOimo GUARD PIPE FROM CENTER OF dRtAlv (MM). 
P=A3S0LUTt PRESSURE ON GUARD PlPt < Ml -a ) . 

SPaC£=sPACE BETWEEN PlPtb (MM). 

AMH1F.NI PKESSUHt IS 0.101 MPa (14. 7 0 MSI). 
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TAHLfc 6 

STEADY FLOW PRESSURES ALONG INSIDE OF GUARD PIPE DUE TO STEAM JET 

FROM BREAK IN PROCESS PIPE 

Z=U1STaNCE ALUNb (iUARD PIPE FROM CENTER OF bRt AK (MM> # 
P=AdSOLUTE PRESSURE Oi\ GUARD PIPE (MPA). 

PO=PPESSURE IN PROCESS PIPE (MPA). 

PROCESS PIPE DiAMtTER = 812.d MM <3? IN). 

GUARD PIPE OlAMtTtR = 1086. MM <42.75 IN). 

AMdiENT pressure = o.ioi mpa <14.70 psD. 


PRO'I.NU. 

7 Hr 2 

7 883 

7884 

bRt AK TYPE 

CIRC 
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PU 

6.MRS 
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8.618 
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Figure 17 shows the pressure distributions (indicating the presence of 
the stationary shock) along the plane of symmetry passing through the center 
of the longitudinal break for selected process/guard pipe separation distances. 
The stationary shock (in the steam) shows the characteristic smoothing-out of 
the discontinuity that is inherent in the artificial-viscosity method for 
handling shock discontinuities numerically. 

Some details cf the steady flow (for longitudinal-break configuration) 

In Runs 7 and 11 are presented in Figures 18 and 19, respectively. These 
runs differ only in the spacing between the pipes (4 and 6 inches, respectively) 
The solid curves are isobars. There is a region of backflow towards the break 
in both figures. The entire flow field in Figure 18 is composed of steam. 

Figure 19 shows a trapped streamer of air, which conveniently delineates the 
boundary of the steam jet. 
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658. 



FIGURE 18. FLOW PATTERN FOR PROBLEM NO. 7874 (RUN 7) 
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FIGURE 19. FLOW PATTERN FOR PROBLEM NO. 7876 (RUN 11) 
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CHAPTER 5 
CONCLUSIONS 


The results in Figures 12 - 15 indicate that the steady-flow maximum 
pressure loads on the guard pipe inside surface produced by ruptures 
(longitudinal and circumferential) of the steam process pipe are less than 
the Initial process pipe steam pressure (as Bernollis 1 law would indicate). 

As the process/guard pipe separation distance approaches zero the pressure 
loading approaches the initial process pipe steam pressure. The guard pipe 
pressure loads fall off rapidly with distance from the process pipe break 
site along the guard pipe surface. Transient airshocks with high peak 
pressures (and very short durations) occur before the flow becomes steady; 
however, these airshocks do not transmit sufficient momentum to the guard 
pipe wall to affect the guard pipe structural response appreciably. A guard 
pipe designed to withstand the process pipe pressure can withstand a steam 
jet from a break in the process pipe. The effect of whipping of the process 
pipe is not considered in this report. 

Only simple size scaling can be accurately applied to these results 
(i.e., multiplying all dimensions and time by the same factor). The maximum 
pressure on the guard pipe depends strongly on the ratio of distance to guard 
pipe divided by break width, and weakly on the ratio cf break length (i.e., 
pipe thickness) to break width and on the equation of state for steam. 
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APPENDIX A 

OPENING TIME OF CIRCUMFERENTIAL BREAK IN PROCESS PIPE 


Assume that the break occurs instantaneously and calculate how long it 
takes for the ends to pull apart to the desired break size. The mechanism 
for pulling the ends apart is a relief wave that propagates along the pipe 
and drops the meridional stress from its loaded condition to zero stress. 
The meridional stress is 


0 ' = 36 - 20 MP a s 5250 psi 

where p = pressure in pipe (=8.274 MPa = 1200 psia) 
R = radius of pipe (=0.445m = 17.5 in) 
t 3 thickness of pipe (=50.0 mm = 2 in) 


The speed of a longitudinal wave in an infinite plate is (Reference A1 - p. 81) 


L fp ( 1 - 


T = 5.24 x 10 3 m/s 
v Z ) - 


where C = Young's modulus (=2.0 x 10 5 MPa for steel (Reference A2 - p. 2-68)) 

v = Poisson's ratio (=0.3) 

p = density (=8 x 10 3 kg/m 3 ) 


The speed of the pipe (particle velocity) in the stress-relieved region is 


^Ko'lsky, H. > 

A2 Gray, D. E. 
1972. 


v 



0.948 m/s 


Stress Waves in Solids , Dover, New York, 1963. 

(ed.), American Institute of Physics Handbook , 3rd ed., McGraw-Hill, 


A-l 
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To open a 50.8 min (2 in) break, each end has to move 25.4 mm (1 in). The 
time this takes is 

t (25.4 mm s 1 in) • (25.4 mm) / (0.948 x 10 3 mm/s) ■ 26.8 ms 


In this time, the elastic relief wave has moved down the pipe a distance of 
(0.0268 s) x (5.24 x 10 5 m/s) ■ 140 m. Anything within 70 m that inhibits 
longitudinal motion of the pipe would send back a wave that would slow the 
opening of the gap before it reached a total size of 50.8 mm (2 in). 


A-2 
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